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SIMULATIONS FOR TERRESTRIAL PLANETS FORMATION 
Niu ZHANG 1 ' 2 , Jianghui JI 2 ' 3 
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ABSTRACT 

In this paper, we investigate the formation of terrestrial planets in the late stage of planetary 
formation using two-planet model. At that time, the protostar has formed for about 3 Myr and the 
gas disk has dissipated. In the model, the perturbations from Jupiter and Saturn are considered. We 
also consider variations of the mass of outer planet, and the initial eccentricities and inclinations of 
embryos and planetesimals. Our results show that, terrestrial planets are formed in 50 Myr, and the 
accretion rate is about 60% — 80%. In each simulation, 3-4 terrestrial planets are formed inside 
"Jupiter" with masses of 0.15 — 3.6M®. In the 0.5 - 4 AU, when the eccentricities of planetesimals 
are excited, planetesimals are able to accrete material from wide radial direction. The plenty of water 
material of the terrestrial planet in the Habitable Zone may be transferred from the farther places by 
this mechanism. Accretion could also happen a few times between two major planets only if the outer 
planet has a moderate mass and the small terrestrial planet could survive at some resonances over 
time scale of 10 8 yr. In one of our simulations, com-mensurability of the orbital periods of planets is 
very common. Moreover, a librating-circulating 3:2 configuration of mean motion resonance is found. 
Subject headings: astrophysics-exoplanet-planetary formation-n-body simulation 



1. INTRODUCTION 

Since the discovery of the first extrasolar planet around 
Solar- Type star, the detection of the extrasolar planets 
develops rapidly. To date, more than 300 planets are 
found orbiting their center stars beyond our solar sys- 
tem, including 35 multiple planetary systems. Recently, 
scientists have found evidences of methane and carbon 
dioxide in the atmosphere of a Hot- Jupiter (HD 189733 
b) (http://planetquest.jpl.nasa.gov). Of ~ 300 known 
extrasolar planets, the minimum mass is generally several 
Jupiter masses. There are also several terrestrial planets 
(Super-Earth), but their orbital characteristic is unsuit- 
able for the formation or development of life. Along with 
the development of survey techniques and incoming high 
definition space missions, people will definitely discover 
more and more Earth-like planets in the extrasolar plan- 
etary systems. The research of formation and evolution 
of the terrestrial planet now becomes important topics 
in astrophysics, astrobiology, astroche mistry and so on . 

Planet formation has certain order (jZhou et alJ l2~"005). 
and Jupiter-like planets at greater distance are formed 
faster than those near the Sun. It is generally believed 
that the planet formation may experience the following 
stages: The grains condensed in the initial stage grow 
to km-sized planetesimals in the early stage, and then, 
in the middle stage, Moon-to-Mars sized embryos are 
formed by accretion of the planetesimals. The size of 
embryos correlates with the feeding zone of the plan- 
etesimals. According to the formula of Hill radius: 
R H = r(m/3M Q ) 1 / 3 (where r,m are the heliocentric dis- 
tance and the mass of planetesimal), the distant plan- 
etesimals have wider feeding zones, so the formed em- 
bryos are larger. When the embryos grow to the core of 
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mass (~ 10M$ ), runaway accretion may take place ac- 
cordingly. With more atmosphere accreted, the embryos 
contract, growing ever denser and more mass ive, even- 
tually collapse to for m giant Jovian planets (|Hu fc Xul 
2008: llda fc Lir]l200l . However, at the same period, the 
inner planetesimals accrete in respective accretion scope, 
and then the embryos of terrestrial planet (namely kind 
of the terrestrial planet core matter) are formed. At the 
end of the third stage, it is around that the protostar has 
formed for about 3 Myr, the gas disk has dissipated. A 
few larger bodies with low e and i are in crowds of plan- 
etesimals with certain eccentricities e, and inclinations i. 
In the late stage, the terrestrial planetary embryos are 
excited to high eccentricity orbit by gravitational per- 
turbation. Then, the orbital crossing makes the plan- 
ets accreting material in the broader radial area. Solid 
residue is either scattered out of the planetary system or 
accreted by the massive planet. However, it also has the 
possibility of being ca ptured at the resonance position 
of th e major planets (jNagasawa fc Idal l2000t iHu fc Xul 
l2008h . 

Tak ing our Solar System as the background. IChambersI 
( 2001) made a study of terrestrial planet formation in the 
late stage by numerical simulations. He set 150 — 160 
Moon-to-Mars size planetary embryos in the area of 
0.3 — 2.0 AU, include gravitational perturbations from 
Jupiter and Saturn. He also examined two initial mass 
distributions: approximately uniform masses, and a bi- 
modal mass distribution. The results show that 2 — 4 
planets are formed in 50 Myr, and finally survive over 
200 Myr timescale. The s pace distri b ution and concen- 
tration (see Section 4 in IChambersI (|2001f )) of planets 
formed in the simulations are similar to our solar system. 
However, the planets produced by the simulations usu- 
ally have eccentric orbits wit h higher eccentricities e, and 
inclinations i than Venus. iRavmond. Quinn fc Lunind 
(|2004 I2006D also studied the formation of terrestrial 
planets. In the simulations, they took into account 
Jupiter's gravitational perturbation, wider distribution 
of material (0.5 — 4.5 AU) and higher resolution. The 
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results confirm a leading hypothesis for the origin of 
Earth's water: they may come from the material in the 
outer area by impacts in the late stage o f plan et forma- 
tion. iRavmond. Mandell. fc Sigurdssonl (|2006l ) explored 
the planet formation under planetary migration of the 
giant. In the simulations, super Hot Earth form inte- 
rior to the migrating giant planet, and water-rich, Earth- 
size terrestrial planet are present in the Habitable Zone 
(0.8 — 1.5 AU) and can survive over 10 s yr timescale. 

In our model, Solar System is taken as the background. 
But several changes are worth noting : 1) we use two- 
planet (Jupiter and Saturn, see Fig.l) model. 2) in the 
model, Jupiter and Saturn are supposed to be formed at 
the beginning of the simulation, with two swarms of plan- 
etesimals distributed among 0.5 — 4.2 AU and 6.2 — 9.6 
AU respectively. 3) The initial eccentricities and inclina- 
tions of planetesimals are considered. 4) The variations 
of the mass of Saturn are examined. 5) The exchange of 
material in the radial direction is also studied by the pa- 
rameter of water mass fraction. 6) We perform the simu- 
lations over longer timescale (400 Myr) in order to check 
the stability and the dynamical structure evolution of 
the system. Our results show that the terrestrial planets 
produced interior to Jupiter have higher mass accretion 
rate, and share the similar architecture as the Solar Sys- 
tem. However, the structure beyond Jupiter correlates 
with the initial mass of Saturn. Almost each simulation 
has a water-rich terrestrial planet in the Habitable Zone 
(0.8-1.5 AU). 

In Section 2, the initial conditions, algorithm and in- 
tegration procedure are described in detail. Section 3 
presents the main results. We conclude the outcomes in 
Section 4. 

2. MODEL 
2.1. Initial conditions 

Generally, the time scale for formation of Jupiter-like 
planet is less than 10 Myr (Briccno 2001). Nevertheless, 
as we know, the formation scenario of planet embryos 
is related to their heliocentric distances and the initial 
mass of the star nebular. If we consider the model of 
1.5 MMSN (minimum mass solar nebular), the upper 
bound of the time scale for Jupiter-like planet formation 
correspond s to the time scale fo r the embryo formation 
at 2.5 AU (jKokubo fe Idall2002f) . which is just at 3 : 1 
resonance location of Jupiter. In the region 2.5 — 4.2 
AU, embryos will be cleared off by strong gravitational 
perturbation arising from Jupiter. There should be some 
much smaller solid residue among Jupiter and Saturn, 
even though the 'clearing effect' may throw out most of 
the material in this area. That's why we set embryos 
only in the region 0.5 — 2.5 AU and planetesimals in the 
0.5 - 4.2 AU and 6.2 - 9.6 AU. 

We adopt the surface densit y profile as follows 
(jRavmond. Quinn fc Luninel [20041) : 
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In (1), S s „oM, = 10 g/cm 2 is the surface density 
at snowline, the snowline is at 2.5 AU with Si = 
4 g/cm 2 . As mentioned earlier, the mass of planetary em- 
bryos is proportional to the width of the feeding zone, 
which is associated with Hill Radius, Rh, so the mass of 



an embryo increases as 

M em bryo OC r^(r)R H (2) 

The embryos among 0.5 — 2.5 AU are spaced by A (A 
varying randomly between 2 and 5) mutual Hill Radii, 
Rh,™, , which is defined as 

i^ = ( ^ )(! W )1/3 ( 3 ) 

where a\.2 and TOi,2 are the semi-major axes and masses 
of the embryos respectively. Replacing Rh in (2) with 
RH.m 7 and substituting (1) in (2), then, we get relations 
between the mass of embryos and the parameter A as 

M embryo oc rZ(r)R H , m oc r ^ 4 A 3 ^ 2 (4) 
As shown in Fig. 1, the initial planetesimals are spread 
over 0.5 — 9.6 AU (excluding 3 Hill Radii around Jupiter); 
the distribution of them should meet (1). Here, we 
equally set the masses of planetesimals inside and outside 
Jupiter respectively as shown in (5). Consequently, the 
number distribution of the planetesimals is only needed 
to satisfy TV oc r -1 / 2 . Additionally, we keep the total 
number of planetesimals and embryos inside Jupiter, and 
the number of planetesimals outside Jupiter both equal 
to 200. 

^embryo H~ ^planetesimal, r<,rj u piter 

-^planetesimal . r~>rj uv n er — 200, 
-^-embryo ~r" ^^-planetesimal, r<,rj u pi 
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M planetesimal, r>r,j up iter 

= 7.5M e /200. 

The water mass fraction of th e bodies is same as 
IRavmond. Quinn fc Luninel (|2004[ ), i.e., the planetesi- 
mals beyond 2.5 AU have 5% water material by mass, 
those between 2 — 2.5 AU have 0.1% water material 
by mass, and the others have 0.001% water material 
by mass. The eccentricities and inclinations vary in 
(0 - 0.02) and (0 - 0.05°), respectively. The mass of Sat- 
urn in simulations la/lb, 2a/2b and 3a/3b are O.5M0, 
5M0, 5OM0 respectively. Each simulation is carried 
out twice with a) considering, b) not considering self- 
gravitation of planetesimals among Jupiter and Saturn. 

2.2. Algorithm 

In regular coordinat e system, the motion equa tions of 
an n-body system are (Mu rray fc Dermottlll999D 



(6) 



where the index i = 1, 2, ■ • • , n denotes the body i , and Xi 

, Pi are the generalized coordinate and momentum of the 

n 2 

body i, respectively. The Hamiltonian, H = *Srr 
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G^wij Y — , is the sum of the kinetic and potential 

energy for the system. From (6), we know that the rate 
of any quantity, q , can be conveniently expressed in the 
following form, 



dq v-i r dq dxj , dq dpi \ 

dt ~ l^i Kdxi dt dpi dt I 
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f^^-dxi dpi dpi dxi' 



(7) 
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If we define an operator F = E(£rM ~ wMJ 

i—l 

((Chambers! I1999D . then we can rewrite (7) as = Fq. 
Integral the differential equation over time t\ — t 2 {t 2 > 
ti), and so we get 

q2 = e^ F qi , (8) 

where qi and q 2 are the values of q corresponding to the 
time t\ and t 2 respectively. If we define h = t% — t\ (h 
actually is the internal time step), then expand (8) at 
zero, we have 

h 2 F 2 

q 2 = (l + hF+ — + •••)«!, (9) 

The symplectic integrator is to divide H into pieces, each 
piece could be individually solved, and then they approx- 
imate the solution of the problem via applying the solu- 
tions once a time. For example, we split the Hamiltonian 
H = H\ + H 2 , and hence have the operators Fi,F 2 . It 
is easy to obtain q 2 = e h ( Fl+F2 ' > qi from (8). By expand- 
ing the exponential (attn. F\F 2 ^ F 2 Fi), ignoring the 
second- and higher-order small quantities of h , then we 
have 

e h(F 1 +F 2 ) _ e hFi e hF 2 _|_ h 2 (F 2 F 1 -F 1 F 2 ) 

= 1 + h( F, + F 2 ) + kl£l±^±£il + . . . 

(10) 

We can get a second-order integrator q 2 = 
e hF 2 /2 e hF le hF 2 /2q^ applying a small equivalence 
transformation. 

The key point for symplectic algorithm is how to split 
Hamiltonian H into pieces. Considering a dynamical sys- 
tem composed of N bodies orbiting a massive central 
body, we can split the Ha miltonian H into the primary 
and the secondary parts. [Chambers] (|1999f ) proposed a 
hybrid symplectic algorithm, in which the Hamiltonian 
is divided into the following parts: 

( N 2 

tt ^ / Pi _ Gm G mj \ 

■"1 L^\2m i r iQ /' 

i—l 

N N 

< H 2 =-G£ E ^ (11) 

i=l j=i+l 13 
N 

v u i=i 

where Hi is the unperturbed Keplerian motion of the 
N smaller bodies, H 2 is the total interaction potential 
energy of the N smaller bodies, and H3 is the kinetic 
energy of the center body (Note That: TV has different 
meaning from n above, N refers to the numbers of bod- 
ies excluding the central body). The term of 'hybrid' 
means that, for the convenience of calculation, heliocen- 
tric coordinate s and barycentri c velocities are used while 
solvin g (11) (|Chambersl Il99l |D uncan. Levison &: Led 
1998). From (10), we can infer that all of the higher 
terms depend on both F\ and F\. If F 2 ~ eF\ (e = 
rrii /m a ), and therefore, the second-order integrator 

q 2 = e hF 2 /2 e hF 3 /2 e hF le hF 3 /2 e hF 2 /2 qi ig correct tQ Q^S} 

only when the different parts of Hamiltonian meet the 
conditions Hi ^> H 2 ,Hi ^> H 3 . However, when close 
encounter occurs between two bodies, the distance be- 
tween them Tij approa ches zero, h ence the H\ ^> H 2 
can not be satisfied. Chambers (1999) introduced a 



changeover function K{rij) to translate part of H 2 asso- 
ciated with the close encounter t o H\ , and then integrat e 
it using Bulirsch-Stoer method (Stoer & Bulirsch 1980). 
The modified H\ , H 2 are present as 

( N 2 n 

FT — ( Pi Gm, e mj x 
1 ~ ^ [ 2m,i r i@ I 

i—l 

N N 

-G£ E ^[1-K( nj )], (12) 

N N 

H2 =-ge E ^K( rij ). 

i=lj=t+l 3 

K(rij) tends to zero when rg is small , while tending to 

one when ry is large (|Chamberslll99l . 

We use the hybrid symplectic integrator (|Chamberd 
1999) in MERCURY package to integrate all the simula- 
tions. We take into account that collision and coalescence 
will occur, when the minimum distance between any of 
the two objects is equal to or less than the sum of their 
physical radii. While they were separated by not more 
than 3 Hill radii, we consider close encounters will take 
place. When the distance from the central star is more 
than 100 AU, these bodies are removed, because they are 
so far from the central star that they play an insignifi- 
cant role of the interaction. In addition, we adopt 6 days 
as the length of time step, which is a twentieth period 
of the innermost body at 0.5 AU. The 6 simulations are 
carried out over 400 Myr time scale. At the end of the 
intergration, the changes of energy and angular momenta 
are 10~ 3 and 10 -11 respectively. The 6 simulations are 
performed on a workstation composed of 12 CPUs with 
1.2 GHz, and each costs roughly 45 days. 

3. RESULTS 

All of 6 simulations exhibit some classical processes 
on planet formation. Firstly, we will analyse simulation 
2a/2b to discuss the physical processes which can apply 
to every simulation. Next, we will make a statistical anal- 
ysis in order to find out that how the planet formation 
may rely on different physical factors. 

3.1. Simulation 2a/2b 

At the end of the calculation, 3 — 4 terrestrial plan- 
ets are formed in 2a/2b. Table 1 shows the properties 
of the terrestrial planets from simulations 2a and 2b. 
We label the planets as b, c, d, and so on, according 
to the heliocentric distance (hereinafter). The masses of 
the terrestrial planets range from several Mars masses to 
several Earth masses. All of them are water-rich, except 
the planet e in simulation 2b. Some parameters of a cer- 
tain planet are comparable with the terrestrial planets 
in solar system. For example, the orbital eccentricity of 
planet b in simulation 2b is 0.0309, which is very close 
to that of Earth. 

Fig. 2 is a snapshot of simulation 2a. At 0.1 Myr, it 
is clear that the planetesimals are excited at the 3 : 2 
(3.97 AU),2 : 1 (3.28 AU) and 3 : 1 (2.5 AU) resonance 
locations with Jupiter, and this is similar to the Kirk- 
wood gaps of the asteroidal belt in solar system. For 
about 1 Myr, planetesimals and embryos are deeply in- 
termixed, most of the bodies have large eccentricities. 
Collisions and accretions emerge among planetesimals 
and embryos. This process continues until about 50 Myr, 
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the planetary embryos are mostly formed, and then dy- 
namical evolution is start. The forma tion time scale i n 
our work is in accordance with that of (jlda &: Linl l2004) . 
Finally, inside Jupiter, 3 terrestrial planets are formed 
with masses of 0.15 — 3.6M®. However, at the outer 
region, planetesimals are continuously scattered out of 
the system at 0.1 Myr. For about 10 Myr, there are 
no survivals except at some resonances with the giant 
planet. As shown in Fig. 2, there is a small body at the 
1 : 2 resonance with Jupiter. Due to the planetesimals' 
scattering, Jupiter (Saturn) migrates inward (outward) 
0.13 AU (1.19 AU) toward the sun respectively Such 
kin d of migration agrees with the work of Fernandez et 
al. (|Fernandez Sz. Idll984[ ) Hence, the 2 : 5 mean motion 
resonance is destroyed, then the ratios of periods between 
Jupiter and Saturn degenerate to 1 : 3. Therefore, the 
ratio of periods for Jupiter, small body and Saturn is 
approximate to 1 : 2 : 3. 

Fig. 3 is a snapshot of simulation 2b. In comparison 
with Fig. 2, it is apparent that planetesimals are excited 
more quickly at the 3 : 2 (3.97 AU), 2 : 1 (3.28 AU) and 

3 : 1 (2.5 AU) resonance location with Jupiter. The sev- 
eral characteristic time scales are the same as simulation 
2a for the bodies within Jupiter. 4 planets are formed 
in simulation 2b, the changes of position of Jupiter and 
Sat-urn are about the same as simulation 2a. We note 
simulations 2a and 2b have the same initial conditions, 
the only difference between them is whether we consider 
the self- gravitation among the outer planetesimals. The 
results of simulation 2b are shown to be a consequence 
of being expected but not surprising. There is a little 
stack planetesimals survival over 400 Myr among 7 — 8 
AU, located in the area of 2 : 3 (6.63 AU) and 1 : 2 (8.03 
AU) resonances with Jupiter. 

Planets in 55 enc planetary syste m have similar spa - 
tial distribution to the solar system (|Fischer et al.l [2008). 
From Table 1 and Fig. 3, it is not difficult to see that 

4 terrestrial planets formed in simulation 2b move on 
the nearly-circular orbit. Mars is ever regarded as a sur- 
vivor of an original planetary embryo, according to its 
unique chemical and isotopic characteristics. As a mat- 
ter of fact, the planet e in simulation 2b is a survivor 
of the initial planetary embryos. Planet e in simulation 
2b does not accrete anything over 400 Myr integration 
time. Comparing the semi-major axis of Mars with that 
of planet e, we notice that the planetesimals of simula- 
tion 2b are located in the asteroidal belt. Furthermore, 
the ratio of periods of the planet e and Jupiter is nearly 
1:2. The total mass of the main-belt i n solar system is 
about 5 x 10- 5 M®, being 0.1% - 0.12% dHu fc Xull20f)l 
of the initial solid material. If the assumed planet e in 
simulation 2b would break into thousands of fragments, 
they may undergo re-accretion or ejection by the pertur- 
bation of Jupiter over secular evolution. If it is similar to 
the same ratio of mass of the belt in solar system, then 
the leftovers of the solid materials almost bear a total 
mass of 6.3 x 10~ 5 — 7.56 x 10~ 5 M^ . In this sense, an 
asteroidal belt is very likely to form in the system quite 
similar to that of our solar system. 

Fig. 4a is the mass curve of terrestrial planets for sim- 
ulation 2b. We can find that the accrete velocity is not 
uniform. At 10 Myr, planets reach half mass of their 
final mass, and then, the accretion velocity slows down, 
because the planetesimals are only a quarter left. Until 



about 50 Myr, the terrestrial planets are formed. The 
accretion rate and mass concentration (the mass rate of 
the largest terrestrial planet and the total formed ob- 
jects) are 73% and 43%, respectively. The corresponding 
parameters in simulation 2a are 60% and 81% respec- 
tively. However, in the area outside Jupiter, 81% initial 
material is scattered out of the system. 

Planet embryos are formed from feeding zones where 
the planetesimals are located. A feeding zone has unique 
chemical and isotopic characteristics. It is helpful to 
study the trace of the planetesimals to understand the 
composition of terrestrial planet, vice versa, for example, 
if we can investigate the origin and formation process by 
revealing the chemical or isotopic characteristics of the 
moon. In Fig. 4b is shown the trace of survivals. We 
note that all the materials accreted by terrestrial planets 
come from the inner swarm of planetesimals or embryos. 
Here Jupiter is like a wall, which separates the inner and 
outer planetesimals from exchanging materials. Once 
again, Fig. 4b verifies that Jupiter may protect the inner 
terrestrial plane ts from colliding with the outer bodies 
(jWetheril]lll990h . We set a water mass fraction on each 
body, it is easy to work out how much water-material 
of the finally terrestrial planet bears. Take planet c in 
simulation 2b for example, the water material is approx- 
imately 1.1 x I0 22 kg , about 8 times than Earth. Fig. 4b 
can also verify that terrestrial planet accrete material in 
broad radial direction. 

3.2. Statistical analysis 

The production efficiency of the terrestrial planet in 
our model is high, and the accretion rate inside Jupiter 
is 60% — 80% in the simulations. 3 — 4 terrestrial planets 
formed in 50 Myr. 5 of 6 simulations have a terrestrial 
planet in the Habitable Zone (0.8 — 1.5 AU) (see Fig. 
5). The planetary systems are formed to have nearly cir- 
cular orbit and coplanarity, similar to the solar system 
(see Table 2). We suppose that the above characteristics 
are correlated with the initial small eccentricities and in- 
clinations. Such adoption could generate more close en- 
counters or collisions in the early several Myr, which may 
increase viscosity of the system and then make the orbits 
more nested on circular orbit on the orbital plane. The 
concentration in Table 2 means the ratio of maximum 
terrestrial planet formed in the simulation and the total 
terrestrial planets mass. It represents different capability 
on accretion, and is not associated with self-gravitation. 
The average value of this parameter is similar to the solar 
system. Considering the self-gravitation of planetesimals 
among Jupiter and Saturn, the system has a better vis- 
cosity, so that the planetesimals will be excited slower. 
The consideration of self-gravitation may not change the 
formation time scale of terrestrial planets, but will af- 
fect the initial accretion speed and the eventual accre- 
tion rate. In Table 2, the simulations lb, 2b, 3b have 
a bit higher accretion rate. When the self-gravitation is 
not considered, the planetesimals may be excited quickly. 
The accretion has a faster speed at the early several Myr, 
so this can promote the accretion rate. 

From Fig. 4b, we have to be aware of that Saturn 
accretes a few planetesimals, this is uncommon in sim- 
ulations la, lb, 3a, 3b. And Fig. 5 shows the finally 
structure of the simulations, it is clear that different Sat- 
urn mass will affect the outer structure of the system 
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beyond Jupiter. In simulation 3a (3b), the Saturn mass 
is 50M® . Now it is large enough to clear the area among 
Jupiter and Saturn. In simulation la (lb), Saturn's mass 
is O.5M0 , more or less equal to the embryos' mass. In 
this case, it is too small to clear off any planetesimal 
amongst the region of Jupiter and Saturn. Therefore, we 
draw the conclusion that only the Saturn's mass is close 
to be 5 — 10M® , then accretion may happen. Saturn 
and Jupiter in our solar system may form in same stage. 
If there exist embryos of IOMq outside Saturn, the giant 
Jupiter-mass planets may form. 

The scattering of planetesimals could cause the migra- 
tion of planets, for example, Jupiter migrated from 5.2 
AU to 5.06 AU while Saturn traveled from 9.6 AU to 
10.71 AU. There are plenty of ratios of semi-major axis 
of the survival planets nearly 2:1. In this case, the 
planet is easily to be captu r ed on 2 : 1 mean motion res- 
onance ([Lee fe Peald l200l Heel [200l IZhou et al.1 [20051: 
iFischer et al.ll2008f ). There are still some ratios of peri- 
ods between the survival bodies close to a simple ratio of 
integers (see Fig. 5) . In the very long process of dynam- 
ical evolution after planetary formation, the planets also 
have the possibility of been captured onto resonant or- 
bit. For example, the orbits show the 2 : 3 mean motion 
resonance be-tween Jupiter and outer small body in sim- 
ulation lb, and the 3 : 1 resonance between Jupiter and 
the planet d in simulation 3a and so on. Many researchers 
have stud ied the resonance and stabili t y of the planetary 
syste ms (|Ji et al.1 l2003t IZhou fc Sur] (20031 : IZhou et al.l 
|2004| ). As shown in Fig. 6, in simulation 2a, two planets 
are on crossing orbits. When a close encounter occurs, 
3 : 2 mean motion resonance is formed, with resonance 
angle 6\ = 2Ai — 3A2 + vji (where Ai t 2 are the mean lon- 
gitude and the longitudes of periapse, the footnotes 1, 2 
means the inner and outer planets respectively.) librating 
around 180° , while O2 — 2Ai — 3A2 +td\ circulating, and 
ei shows large oscillations. Such 'librating-circulating' 
configuration is similar to the configuration of 2 : 1 res- 
onance in HD 73526 planetary syste m. There have been 
several hypotheses about its origin (Tinney et al. 120061 : 
ISandor fc Klevl [20061 : ISandor. Kiev fc Klagyivikl 120071 ). 
However, it still needs further study in the future. 

4. CONCLUSIONS 

We simulate the terrestrial planets formation by using 
two-planet model. In the simulation, the variations of 
the mass of outer planet, the initial eccentricities and 
inclinations of embryos and planetesimals are also con- 
sidered. The results show that, during the terrestrial 
planets formation, planets can accrete material from dif- 
ferent regions inside Jupiter. Among 0.5 — 4.2 AU, the 
accretion rate of terrestrial planet is 60% — 80%, i.e., 
about 20% — 40% initial mass is removed during the 
progress. The planetesimals will improve the efficiency 
of accretion rate for certain initial eccentricities and in- 
clinations, and this also makes the newly-born terrestrial 
planets have lower orbital eccentricities. It is maybe a 
common phenomenon in the planet formation that the 



Briceno, C. et al., 2001, Science, 291, 93 
Chambers, J. E. 1999, MNRAS, 304, 793 
Chambers, J. E. 2001, Icarus, 152, 205 

Duncan, M., Levison, H. F., & Lee, M. H. 1998, AJ, 116, 2067 



water-rich terrestrial planet is formed in the Habitable 
Zone. The structure, which is similar to that of solar 
system, may explain the results of disintegration of a ter- 
restrial planet. Most of the planetesimals among Jupiter 
and Saturn are scattered out of the pl anetary systems, 
and t his migration caused by scattering ([Fernandez fc Ipl 
1984) or long-term orbital evolution can make planets 
capture at some mean motion resonance location. Accre- 
tion could also happen a few times between two planets 
if the outer planet has a moderate mass, and the small 
terrestrial planet could survive at some resonances over 
10 8 yr time scale. Structurally, Saturn has little effect 
on the architecture inside Jupiter, owing to its protec- 
tion. However, obviously, a different Saturn mass could 
play a vital role of the structure outer Jupiter. Jupiter 
and Saturn in the solar system may form over the same 
period. 

In our simulations, neither terrestrial planets are 
formed within 0.1 AU, nor planetesimals or embryos are 
left. However, a lot of exoplanets with orbital semi- 
major less than 0.1 AU are observed, and several Super- 
Earths are discovered. It is usually believed that they 
were formed far from the center star and then migrated 
into current location ([Raymond. Mandell. fc Sigurdssonl 
2006). We do not consider the migration in the sim- 
ulations, which is caused by the interaction between 
the giant plane ts or planetesimals in the gaseous disk 
(jOu et al.ll2007t ). So the simulation in this work can be 
applied to the case of the dissipation of gas disk, in the 
late stage of planet formation. In the future study, we 
will consider the giant planets under inward migration, 
and in such circumstances short-period terrestrial planet 
could be produced. To date, terrestrial planets are not 
detected in the observations, due to the reasons of the 
selection effect of detection methods and the low res- 
olution precision. Both Doppler velocities and transit 
method are sensitive to the objects moving in smaller or- 
bits. The current research of extrasolar terrestrial plan- 
ets has greatly contributed to the origin and evolution of 
our own solar system. Kepler has been launched success- 
fully on March 6, 2009, whose main scientific objective is 
detecting the Earth- like terrestrial planets. Along with 
high accuracy incoming space projects, it is predictable 
that more and more extrasolar planetary systems with 
similar structure to the solar system will be discovered. 
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TABLE 1 

Properties of terrestrial planets from simulations 2a and 2b. 



Planet 


a (AU) 


e 


i(deg) 


m(m©) 


water mass 


2a b 


0.6174 


0.1098 


7.452 


3.6421 


0.0316% 


c 


1.7796 


0.1935 


33.890 


0.1528 


0.1040% 


d 


2.3304 


0.3291 


9.393 


0.6925 


0.5218% 


2b b 


0.5274 


0.0309 


4.041 


2.3527 


0.1539% 


c 


1.0451 


0.1080 


4.915 


1.3880 


0.1316% 


d 


1.4783 


0.0520 


5.189 


1.6696 


0.7795% 


e 


3.1162 


0.2086 


6.421 


0.0630 


0.0010% 



TABLE 2 

Properties of terrestrial planets from different systems 



System 


accretion rate 


n 


ro(me) 


concentration 


e 


in 


la 


73.2518% 


3 


1.8313 


0.4606 


0.1381 


7.6963 


lb 


80.3853% 


3 


2.0096 


0.4262 


0.0937 


1.7790 


2a 


59.8322% 


3 


1.4958 


0.8116 


0.2108 


16.9117 


2b 


72.9779% 


4 


1.3683 


0.4299 


0.0999 


5.1415 


3a 


65.1098% 


3 


1.6277 


0.5337 


0.2063 


5.9153 


3b 


66.9694% 


3 


1.6742 


0.5040 


0.1839 


5.2447 


la-3b 


69.7544% 


3.2 


1.6678 


0.5276 


0.1554 


7.1148 


solar 




4 


0.4943 


0.5058 


0.0764 


3.0624 
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Fig. 1. — Two-planet model and initial conditions. 
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Fig. 2. — Snapshot of simulation 2a with Msatum = 5Mg5. The total mass of embryos is 2.4Mg, the masses of planetesimals inside 
Jupiter are O.O317M0, and those outside Jupiter are O.O375M0. Planetesimals among Jupiter and Saturn were nonself-gravitational (see 
Section 2.1). Note the size of each object is relative, and the value bar is log of water mass fraction, e.g. the wettest body has water mass 
fraction log 10 (5%) = —1.3. 
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Fig. 4. — (a) Mass curve in simulation 2b. 




Fig. 4. — (b) Trace of the objects' in simulation 2b. 
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Fig. 5. — Results of six simulations. The • and X denote formed terrestrial planets and survival planetesimals. Jupiter and Saturn locate 
at about 5 AU and 10 AU, respectively. H zone marked with dotted lines is so called Habitable Zone in 0.8 - 1.5 AU. Some mean motion 
resonance locations with Jupiter are also labeled in the figure. 
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Fig. 6. — Resonant variable 81 of the 3:2 mean motion resonance of two terrestrial planets and the curve of eccentricities and semi-major 
axes in simulation 2a. (6\ = 2Ai — 3A2 + where Ai^ are the mean longitudes and the longitudes of periapse, the footnotes 1, 2 mean 
the inner and outer planets respectively.) 



